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Abstract 

Nanomechanical resonators (NRs) with very high frequency have a great potential for mass sensing with 
unprecedented sensitivity. In this study, we propose a scheme for mass sensing based on the NR capacitively 
coupled to a Cooper-pair box (CPB) driven by two microwave currents. The accreted mass landing on the 
resonator can be measured conveniently by tracking the resonance frequency shifts because of mass changes in 
the signal absorption spectrum. We demonstrate that frequency shifts induced by adsorption often 1587 bp DNA 
molecules can be well resolved in the absorption spectrum. Integration with the CPB enables capacitive readout of 
the mechanical resonance directly on the chip. 



1 Introduction 

Nanoelectromechanical systems (NEMS) offer new pro- 
spects for a variety of important applications ranging 
from semiconductor-based technology to fundamental 
science [1]. In particular, the minuscule masses of 
NEMS resonators, combined with their high frequencies 
and high resonance quality factors, are very appealing 
for mass sensing [2-7]. These NEMS-based mass sensing 
employs tracking the resonance frequency shifts of the 
resonators due to mass changes. The most frequently 
used techniques for measuring the resonance frequency 
are based on optical detection [8]. Though inherently 
simple and highly sensitive, this technique is susceptible 
to temperature fluctuation noise because it usually gen- 
erates heat and heat conduction. On the other hand, it 
has experimentally been demonstrated that capacitive 
detection is less affected to noise than optical detection 
in ambient atmosphere [9]. Capacitive detection is rea- 
lized by connecting NEMS resonator with standard 
microelectronics, such as complementary metal-oxide- 
semiconductor (CMOS) circuitry [10]. Here, we propose 
a scheme for mass sensing based on a coupled nanome- 
chanical resonator (NR) -Cooper-pair box (CPB) system. 

The basic superconducting CPB consists of a low- 
capacitance superconducting electrode weakly linked to 
a superconducting reservoir by a Josephson tunnel 
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junction. Owing to its controllability [11-14], a CPB has 
been proposed to couple to the NR to drive an NR into 
a superposition of spatially separated states and probe 
the decay of the NR [15], to prepare the NR in a Fock 
state and perform a quantum non-demolition measure- 
ment of the Fock state [16], and to cool the NR to its 
ground state [17]. Recently, this coupled CPB-NR sys- 
tem has been realized in experiments [18,19] and the 
resonance frequency shifts of the NR could be moni- 
tored by performing microwave (MW) spectroscopy 
measurement. Based on the above-mentioned achieve- 
ments, in this article, we investigate the signal absorp- 
tion spectrum of the CPB qubit capacitively coupled to 
an NR in the simultaneous presence of a strong control 
MW current and a weak signal MW current. Theoreti- 
cal analysis shows that two sideband peaks appear at the 
signal absorption spectrum, which exactly correspond to 
the resonance frequency of the NR. Therefore, the 
accreted mass landing on the NR can be weighed pre- 
cisely by measuring the frequency shifts because of mass 
changes of the NR in the signal absorption spectrum. 
Similar mass sensing scheme has been proposed recently 
in a hybrid nanocrystal coupled to an NR by our group 
[20], which is based on a theoretical model. However, 
recent experimental achievements in the coupled CPB- 
NR system [18,19] make it possible for our proposed 
mass sensing scheme here to be realized in future. 
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Figure 1 Schematic diagram of an NR capacitively coupled to a CPB. Two MW currents with frequency co c and co s and a direct current / b 
are applied in the MW line to control the magnetic flux & x through the CPB loop. 



2 Model and theory 

In our CPB-NR composite system shown schematically 
in Figure 1, the NR is capacitively coupled to a CPB 
qubit consisting of two Josephson junctions which form 
a SQUID loop. A strong control MW current and a 
weak signal MW current are simultaneously applied in a 
MW line through the CPB to induce the oscillating 
magnetic fields in the Josephson junction SQUID loop 
of the CPB qubit. Besides, a direct current I b is also 
applied to the MW line to control the magnetic flux 
through the SQUID loop and thus the effective Joseph- 
son coupling of the CPB qubit. The Hamiltonian of our 
coupled CPB-NR system reads: 



H = Hqpb + + Hi r 
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where //cpb is the Hamiltonian of the CPB qubit 
described by the pseudospin -1/2 operators o z and o x = o 
+ + <7_. &> q = 4E c (2n g - l)lh is the electrostatic energy dif- 
ference and £j 0 is the maximum Josephson energy. Here, 
E c = e 2 l2Cz is the charging energy with = C b + C g + 
2Cj being the CPB island's total capacitance and n g = 



(CbVb + C g Vg)/(2e) is the dimensionless polarization 
charge (in units of Cooper pairs), where C b and V h are, 
respectively, the capacitance and voltage between the NR 
and the CPB island, C g and V g are, respectively, the gate 
capacitance and voltage of the CPB qubit, and Q is the 
capacitance of each Josephson junction. Displacement 
(by x) of the NR leads to linear modulation of the capaci- 
tance between NR and CPB, C b (x) * C b (0) + (dC b /dx)x, 
which modulates the electrostatic energy of the CPB 
qubit, resulting in the capacitive coupling constant 

k= 4 j£E£^ AXzp , where = C b V h /2e and 



Ax zp = y/h/2mco n is the zero-point uncertainty of the NR 
with effective mass m and resonance frequency co n . The 
coupling between the MW line and the CPB qubit in the 
second term of Equation 2 results from the totally exter- 
nally applied magnetic flux & x (t) = O q (t) + & h through 
the CPB qubit loop of an effective area S with & 0 = h/(2e) 
being the flux quantum. Here, 0 q (t) = (j 0 SI(t)/(2jTr), with 
r being the distance between the MW line and the qubit 
and being the vacuum permeability. & q (t) and & h are 
controlled, respectively, by the MW current 
J(0 = £ c cos(cD c t) + £ s cos(co s t + 8 f ) and the direct current 
Yb in the MW line. For convenience, we assume the phase 
factor d' = 0 because it is not difficult to demonstrate that 
the results of this article are not dependent on the value 
of S\ By adjusting the direct current I h and the MW cur- 
rent I{t) such that <2> b » <2> q (t) and 7i& h /& 0 = n /2, we 

can obtain Ejcos [ ^p ] ^ —E] n0 ^\ In a rotating 
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frame at the control frequency co c , the total Hamiltonian 
can now be written as 

H = —hAa z + hco n a^a + hX{a^ + a)a z + hQ(a+ + a_) 
2 (5) 

+/i£ s {cr + e- i8t + cr-e i8t ), 

where A = co q - co c is the detuning of the qubit reso- 
nance frequency and the control current frequency, S = 
co s - co c is the detuning of the signal current and the 
control current, (A = ju 0 SE 10 / (8r& 0 ) is the effective 'elec- 
tric dipole moment' of the qubit, and Q = jji£ c /h is the 
effective 'Rabi frequency of the control current. 

The dynamics of the coupled CPB-NR system in the 
presence of dissipation and dephasing is described by 
the following master equation [21] 

% = ~[H, P ] + ^rC[*-] + $C[a] + ^-C[a z ], (6) 
at n 2T\ 2 4r^ 

where p is the density matrix of the coupled system, 
Ti is the qubit relaxation time, r<p is the qubit pure 
dephasing time, and y is the decay rate of the NR which 
is given by y - &>«/Q- £[D], describing the incoherent 
decays, is the Lindblad operator for an operator and is 
given by: 



C[D] = IDprf - D^Dp - pD^D. 



(7) 



Using the identity (6) = Tr(Op) for an operator O and 
a density matrix p in Equation 6, we obtain the follow- 
ing Bloch equations for the coupled CPB-NR system: 



d(<7- 
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equations, we first take the semiclassical approach by 
factorizing the NR and CPB qubit degrees of freedom, i. 
e., (Q(7_) = (Q) (cj_), which ignores any entanglement 
between these systems. For simplicity, we define p - ft<j_ 
, k = G z and then we have 



dp 
dt 
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where £ = £ c + £ s e~ m . In order to solve the above 
equations, we make the ansatz (p(t)) = p 0 + pie~' St + p. 
id*, (k(t)) = k 0 + k ie - i3t + k^e iSt , and <Q(f)> = Q 0 + 
lSt + Q.\e lSt [22], Upon substituting these equations into 
Equations 12-14 and upon working to the lowest order 
in £ s but to all orders in £ c , we obtain in the steady 
state: 

H 2 £ s T 2 k 0 2T 1 /T 2 B(<5 0 + 2i)(C + Q 2 C ) + E(B - So) 
Pi = 1 • ( ib > 



AE{B-S 0 ) 



where 



(16) 



A=A C — 4X 0 co 0 k 0 - 8 0 - i, 
B = A c — AXqCjOqIzo + 80 + i, 
C = AX 0 cookor)Qll{A c - 4X 0 co 0 k 0 - i), 
D = 4X 0 co 0 k 0 riQc/i A c ~ ^oco 0 k 0 + i), 
E = 2T l /T 2 A(D + Ql) - 2Ti/T 2 B(C + Q^) 
— AB(Ti/T 2 8 0 + i). 

Here, dimensionless variables co 0 = co r T 2 , /o = 7^2> Q c 
= AT 2) and A c = AT 2 are introduced for convenience 
and the auxiliary function 



Y] 



u>o - iyo8 0 - 8 2 



(17) 



d 2 {Q) d{Q) 2 3 
___ + y __ + , r(Q) = _4 Wr Ao(a z ), 



where A 0 = and T 2 is the qubit dephasing time 
satisfying 



1 1 1 

T 2 2Ti r<p ' 



(11) 



Note that if the pure dephasing rate is neglected, i.e., 
^ = 0, then T 2 = 2T X . In order to solve the above 



The population inversion /c 0 of the CPB is determined 



(10) by 



(feo + 1)[(A C - 4X 0 co 0 k 0 ) 2 + 1] + 4^ c 2 fe 0 - i = 0. (18) 

Pi is a parameter corresponding to the linear suscept- 
ibility x [1 \oj s ) = pxl£ s = (fA 2 T 2 /h)x{co s ), where the 
dimensionless linear susceptibility #(&> 5 ) is given by 

, , m/T^BjSo + 2Q(C + i2?) + £(B - S 0 ) 
X(Ws) = AE{B - S 0 ) feo U9) 
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The real and imaginary parts of %(co s ) characterize, 
respectively, the dispersive and absorptive properties. 

The coupled CPB-NR system has been proposed to 
measure the vibration frequency of the NR by calculat- 
ing the absorption spectrum [23]. On the other hand, 
NRs have widely been used as mass sensors by measur- 
ing the resonant frequency shift because of the added 
mass of the bound particles. The mass sensing principle 
is simple. NRs can be described by harmonic oscillators 
with an effective mass m ef f, a spring constant /c, and a 
mechanical resonance frequency co n = ^k/m e ff. When a 
particle adsorbs to the resonator and significantly 
increases the resonator's effective mass, therefore, the 
mechanical resonance frequency reduces. Mass sensing 
is based on monitoring the frequency shift Acq of co n 
induced by the adsorption to the resonator. The rela- 
tionship between Acq with the deposited mass Am is 
given by 

Am = -^Hh^Aco = K~ l Aco, (20) 

where 1Z = (— 2ra e ff/&> n ) -1 is defined as the mass 
responsivity. However, the measurement techniques are 
rather challenging. For example, electrical measurement 
is unsuitable for mass detections based on very high fre- 
quency NRs because of the generated heat effect [24]. 
For optical detection, as device dimensions are scaled 
far below the detection wavelength, diffraction effects 
become pronounced and will limit the sensitivity of this 
approach [25]. Moreover, in any actual implementation, 
frequency stability of the measuring system as well as 
various noise sources, including thermomechanical noise 
generated by the internal loss mechanisms in the reso- 
nator and Nyquist-Johnson noise from the readout cir- 
cuitry [3,26] will also impose limits to the sensitivity of 
measurement. Here, we can determine the frequency 
shifts with high precision by the MW spectroscopy mea- 
surement based on our coupled CPB-NR system. 

3 Numerical results and discussion 

In what follows, we choose the realistically reasonable 
parameters to demonstrate the validity of mass sensing 
based on the coupled CPB-NR system. Typical para- 
meters of the CPB (charge qubit) are E c lh = 40 GHz 
and Ej 0 /h = 4 GHz such that E c » £j [27] . Experiments 
by many researchers have demonstrated CPB eigenstates 
with excited state lifetime of up to 2 fts and coherence 
times of a superpositions states as long as 0.5 fts, i.e., I\ 
= 2^s, and T 2 = 0.5 {is [13,28,29]. NR with resonance 
frequency co n = 2tt x 133 MHz, quality factor Q = 5000, 
and effective mass m e ff = 73 fg has been used for zepto- 
gram-scale mass sensing [5], Besides, coupling constant 
A between the CPB and NR can be chosen as A = 0.1co n 



= 2tt x 13.3 MHz [16]. We assume S = 1 ^m 2 , r = 10 
^m, and S c = 200/xA[30], therefore, we can obtain [t/h = 
(JoSEj 0 /(8hr(po) * 30 GHzA" 1 and 
fi c = QT 2 = (/x/fr)£ c T 2 = 3. The experiments of our 
proposed mass sensing scheme should be done in situ 
within a cryogenically cooled, ultrahigh vacuum appara- 
tus with base pressure below 10" 10 Torr. 

Firstly, we would show the principle of measuring the 
resonance frequency of the NR in the coupled CPB-NR 
system. Figure 2a illustrates the absorption of the signal 
current as a function of the detuning A s (A s = co s - co q ). 
The absorption (Im(%)) has been normalized with its 
maximum when the control current is resonant with the 
CPB qubit (A c = 0). Mollow triplet, commonly known in 
atomic and some artificial two-level system [31,32], 
appears in the middle part of Figure 2a. However, there 
are also two sharp peaks located exactly at A s = ±co n in 
the sidebands of the absorption spectrum, which corre- 
sponds to the resonant absorption and amplification of 
the vibrational mode of the NR. Our proposed mass 
sensing scheme is just based on these new features in 
the absorption spectrum. An intuitive physical picture 
explaining these peaks can be given in the energy level 
diagram shown in Figure 2b. The Hamiltonian of the 
coupled system without the externally applied current 
can be diagonalized [33,34] in the eigenbasis of 
|±,N±) = | ± ) z 0 gTCV^X^-a) |jsj), with the eigenener- 
gies E ± = ±hl2co q + hco n {N - A 0 ), where the CPB qubit 
states | ±) z are eigenstates of o z with the excited state | +) 
z = \e) and the ground state \-) z = \g), the resonator 
states \N±) are position-displaced Fock states. Transi- 
tions between |-, N_) and |+, (N + 1)+) represent signal 
absorption centered at co c + co n (the rightmost solid line 
in Figure 2b). Besides, transitions between |+, A/+) and 
|-, (N + 1)_) indicate probe amplification (the leftmost 
solid line in Figure 2b) because of a three-photon pro- 
cess, involving simultaneous absorption of two control 
photons and emission of a photon at frequency co c - co n . 
The middle dashed lines in Figure 2a corresponds to the 
transition where the signal frequency is equal to the 
control frequency. Therefore, Figure 2a provides a 
method to measure the resonance frequency of the NR. 
If we first tune the frequency of the control MW cur- 
rent to be resonant with the CPB qubit (co c = co q ) and 
scan the signal frequency across the CPB qubit fre- 
quency, then we can easily obtain the resonance fre- 
quency of the NR from the signal absorption spectrum. 

Next, we illustrate how to measure the mass of the 
particles landing on the NR based on the above discus- 
sions. Unlike traditional mass spectrometers, nanome- 
chanical mass sensors do not require the potentially 
destructive ionization of the test sample, are more sensi- 
tive to large biomolecules, such as proteins and DNA, 
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Figure 2 Scaled absorption spectrum of the signal current as a function of the detuning A s and energy level diagram of the coupled 

system, (a) Scaled absorption spectrum of the signal current as a function of the detuning 4 without landing any masses on the NR. (b) The 
energy level diagram of the CPB coupled to an NR. (c) Signal absorption spectrum as a function of 4 before (black solid line) and after (red 
dashed line) a binding event of ~ 10 functionalized 1587 bp long dsDNA molecules. Frequency shift of 95 kHz can be well resolved in the 
spectrum. Other parameters used are co n = 835 MHz, A 0 = 0.01, 4 = 0, 0 = 5000, 7"-, = 0.25 ys, T 2 = 0.05 fus, and Q c = 3. 



and could eventually be incorporated on a chip [6]. 
Here, we use the functionalized 1587 bp long dsDNA 
molecules with mass m DNA ~ 1659 zg (1 zg = 10 21 g) 
[35], and assume for simplicity that the mass adds uni- 
formly to the mass of the overall NR and changes the 
resonance frequency of the NR by an amount given by 
Equation 19. Figure 2c demonstrates the signal absorp- 
tion as a function of A s before and after a binding event 
of ~ 10 functionalized 1587 bp DNA molecules in the 
vicinity of the resonance frequency of the NR. We can 
see clearly that there is a resonance frequency shift Acq 
= -95 kHz after the adsorption of the DNA molecules 
because of the increased mass of the NR. According to 
Equation 19, we can obtain the mass of the accreted 
DNA molecule: Am = -^Aco = 16590zg, about the 
mass of 10 functionalized 1587 bp long dsDNA mole- 
cules. Therefore, such a coupled CPB-NR system can be 



used to weigh the external accreted mass landing on the 
NR by measuring the frequency shift in the signal 
absorption spectrum when the control current is reso- 
nant with the CPB qubit. Plot of frequency shifts versus 
the number of DNA molecules landing on two different 
masses of NRs. Other parameters used are co n = 835 
MHz, A 0 = 0.01, A c = 0, Q = 5000, T x = 0.25 ^s, T 2 = 
0.05 fts, and O c = 3. Mass responsivity Tl is an impor- 
tant parameter to evaluate the performance of a resona- 
tor for mass sensing. Figure 3 plots the frequency shifts 
as a function of the number of DNA molecules landing 
on the NR for two different kinds of NRs. One is co n = 
2tt x 133 MHz (m e ff = 73 fg), the other is co n = 2tt x 
190 MHz (m eff = 96 fg) [2,3]. The mass responsivities, 
which can be obtained from the slope of the line, are, 
respectively, \1Z\ ~5.72Hz/zg and \TZ\ ~ 6.21 Hz/zg. 
Smaller mass of the nanoresonator enables higher mass 
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Figure 3 Plot of frequency shifts versus the number of DNA molecules landing on two different masses of NRs. Other parameters used 
are co n = 835 MHz, A 0 = 0.01, 4 = 0. 0 = 5000, h = 0.25 jus, T 2 = 0.05 ijs, and O c = 3. 



responsivity. Here, we have assumed that the DNA 
molecules land evenly on the NR and they remain on it. 
In fact, the position on the surface of the resonator 
where the binding takes place is one factor that strongly 
affects the resonance frequency shift. The maximization 
in mass responsivity is obtained if the landing takes 
places at the position where the resonator's vibrational 
amplitude is maximum. For the doubly clamped NR 
used in our model, maximum shift is achieved at the 
center for the fundamental mode of vibration, while the 
minimum shift exists at the clamping points. This statis- 
tical distribution of frequency shifts has been investi- 
gated by building the histogram of event probability 
versus frequency shift for small ensembles of sequential 
single molecule or single nanoparticle adsorption events 
[6,7]. 

In order to demonstrate the novelty of our proposed 
mass sensing scheme, we plot Figure 4 to illustrate how 
the vibration mode of NR and the control current affect 
the spectral features. Figure 4a shows the absorption 
spectrum of the signal field through the CPB system 
without the influence of the NR (coupling off) in the 



absence of the control field (control off), which shows 
the standard resonance absorption profile. However, 
when the coupling turns on, the center of the curve 
shifts from the resonance co s = co q a bit, as shown in 
Figure 4b. This is because of the coupling A 0 between 
the CPB and the NR [16,36]. Figure 4c demonstrates the 
absorption spectrum of the signal field when the control 
field turns on in the absence of the NR (coupling off). 
This is the commonly known Mollow triplet, which 
appears in atomic and some artificial two-level system 
[31,32]. None of the above situations can be used to 
measure the resonance frequency of the NR. However, 
when the coupled CPB-NR system is driven by a strong 
control field and a weak signal field simultaneously, the 
resonance frequency of the NR be measured from the 
absorption spectrum of the signal field, as shown in Fig- 
ure 4d. The spectral linewidth of the two sideband 
peaks that corresponds to the resonance frequency of 
the NR is much narrower than the peak in the center, 
since the damping rate of the NR is much smaller than 
the decay rate of the CPB qubit. Therefore, such a 
coupled CPB-NR system is proposed here to measure 
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Figure 4 Signal current absorption spectrum as a function of the detuning A s considering the effects of NR and the control field 

Other parameters are 4 = 0, Q = 5000, ^ = 0.25 fjs, T 2 = 0.05 fjs, and co n = 835 MHz. 



the resonance frequency of the NR when the control 
field is resonant with the CPB qubit (co c = co q ). By mea- 
suring the frequency shift of the NR before and after 
the adsorption of particles landing on it, we can obtain 
the accreted mass according to Equation 19. 

4 Conclusion 

To conclude, we have demonstrated that the coupled 
NR-CPB system driven by two MW currents can be 
employed as a mass sensor. In this coupled system, the 
CPB serves as an auxiliary system to read out the reso- 
nance frequency of the NR. Therefore, the accreted 
mass landing on the NR can be weighed conveniently 
by measuring the frequency shifts in the signal absorp- 
tion spectrum. In addition, the use of on-chip capacitive 
readout will prove especially advantageous for detection 
in liquid environments of low or arbitrarily varying opti- 
cal transparency, as well as for operation at cryogenic 
temperatures, where maintenance of precise optical 
component alignment becomes difficult. 
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